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^ ' Abstract 

For problems with multiple bodies, the current integral approach 
needs the use of auxiliary potential functions in order to have an in- 
dividual force formula for each body. While the singularity approach, 
based on an extension of the unsteady Lagally theorem, is restricted 
Cn ', to multibody and multivortex flows without bound vortex and vortex 

production. In this paper, we consider multibody and multivortex 
^T . flow and derive force formulas, in both forms of singularity approach 

^f^ . and integral approach but without auxiliary function, that give indi- 

vidual forces of each body for unsteady two dimensional potential flow 
with vortex production on the surface of bodies. A number of prob- 
lems, including Karman vortex street, Wagner problem of impulsively 
starting flow, interaction of two circular cylinders with circulation, 
^.' and interaction of an airfoil with a bound vortex, are used to validate 

the force formulas. 
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1 Introduction 

In the classic Kutta Joukowski theorem, the role of the starting vortex, pro- 
duced during the starting up of flow, is omitted by simply assuming it disap- 
pear in the far flow field. Thus the lift is related to the circulation of bound 
vortex (Batchelor 1967). Under the assumption of steady potential flow, the 
circulation of the bound vortex is determined by the Kutta condition, due to 
which the role of viscosity is implicitly incorporated though explicitly ignored 
(Grighton 1985). The lift predicted by Kutta Joukowski theorem within the 
framework of inviscid flow theory is quite accurate even for real viscous flow, 
provided the flow is steady and unseparated, see Anderson (1984,p.l92) for 
more details. 

For many problems, there may be free vortices, including the starting vor- 
tex, or other bodies close to the body. These problems have been attracting 
great attentions since more than two decades ago, due to their wide applica- 
tions in unsteady flows (Chow&Huang 1982, Lee&Smith 1991, Aref 2007) and 
in multibody flows such as multi-turbine flow (Oterberg, 2010), multi-blade 
flow (Smith & Timoshin, 1996), multi-element airfoil flow (Katz & Plotkin, 
2001), multi-wing aerodynamics as for dragonfly (Hsieh, Kung&Chang 2010), 
and flows in staggered cylinders (Crowdy 2006). Early studies considering 
the interaction of a vortex with a wing can be found in Saffman (1992,pl22). 
The recent studies led to force formulas which can be conveniently classifled 
into integral approaches and singularity approaches. For integral approaches, 
the forces are expressed in terms of the time variation of the integrated vor- 
ticity moment or fluid impulse. For singularity approaches, the forces are 
expressed in terms of the production of the strengths of singularities and the 
speeds of the real and image singularities. In both approaches, there may be 
additional terms, including for instance the added mass effect, if the body is 
subjected to acceleration, rotation and deformation. 

The major differences between various integral approaches lie on the 
choice of the integration domain. The integral approaches of Wu (1981) for 
viscous flow and Saffman (1992) for inviscid flow use integrals deflned for the 
whole space including both the fluid region and the solid region. With the 
help of auxiliary flow potentials, Howe (1995) derived an integral approach 
with a surface integral for the viscous force but with the volume integral only 
deflned for the space occupied by the fluid and with the contribution from 
added mass, pressure by vortices and skin frictions represented by separate 
terms. Such an approach has been subsequently extended to multiple bodies 



by Ragazzo & Tabak (2007) and Chang, Yang & Chu (2008). The advantage 
of the integral approach is that it only requires the knowledge of the velocity 
field and its derivatives. The disadvantage is that the vortices in the far flow 
field, such as starting vortices, must also be taken into account, even when the 
starting vortices are at infinity. Integral approaches for a truncated domain 
are then derived, at the sacrifice of introducing boundary surface integrals 
(Noca et al 1999, Wu et al 2007, Eames et al 2008), see Wu, Lu& Zhuang 
(2007) for more recent advances related to integral approaches and for a 
thorough discussion of the usefulness of the various integral approaches. 

The force formulas by the singularity approach are basically worked out 
through using the complex potential theory and the unsteady Blasius theo- 
rem. The Blasius equation for the general case of unsteady flow and for a 
body in arbitrary motion can be found in Thomson (1968). Streitlien and 
Triantafyllou (1995) derived a force formula for a single Joukowski airfoil sur- 
rounded with point vortices convected freely. Ramodanov (2002) considered 
the motion of a circular cylinder in the presence of N point vortices and the 
forces are expressed in terms of the speeds of both real and image vortices. 
Kanso&Oskouei (2008) extended these formulas to deformable bodies with 
vortex production. The forces are also expressed in terms of the velocities of 
real vortices plus the time variation of an additional integral term represent- 
ing all the effects other than the motion of vortices outside of the cylinder, 
see Shashikanth et al (2002) for a circular cylinder, Borisov et al (2007) for 
a general cylinder, and Michelin& Smith (2010) for problems with vortex 
production. 

For singularity approaches, both the real and image vortices are implicitly 
or explicitly included locally according to their real positions. For the case of 
an airfoil interacting with one outside vortex, Katz and Plotkin (2001, chap- 
ter 6.9) express the force in terms of the induced velocity at the body center. 
The role of real vortex outside of the body is represented by the induced ve- 
locity at the body center. This result was obtained under the lumped vortex 
assumption and extended to the case of multi-airfoil and multi- vortex flow by 
Bai&Wu (2013). The presence of singularities such as sources and doublets 
outside of a body has been also studied in the framework of Lagally theorem 
(Milne-Thomson 1968, Landweber&Milohl980). Wu, Yang &Young (2012) 
extended the Lagally theorem to the case of two dimensional flow with multi- 
body moving in a still fluid in the presence of multiple free vortices, and the 
forces are expressed in terms of the induced velocities or its derivatives at 
the positions of the internal singularities, including sources, doublets and im- 



age vortices. Bound vortex and vortex production are not considered in this 
work. Moreover, validation and application studies are restricted to circular 
cylinders. 

For convenience, the approaches purely based on the velocities of singu- 
larities will be called singularity velocity method. When the induced flow 
velocities at the inner singularities are used to replace the velocities of the 
free singularities, the approach will be called induced velocity method. 

In this paper, we consider unsteady two dimensional potential flow with 
multiple bodies and multiple free vortices. Each body is assumed to have 
an arbitrary shape and vortex production is also considered. Force formulas, 
algebraic and explicit for each body, valid for both discrete and continuously 
distributed singularities, will be derived. This paper will be organized as 
follows. 

In section 2 we first use a momentum approach to relate the lift force and 
induced drag force to the speeds of singularities inside and outside of a single 
body (singularity velocity method) with and without vortex production. We 
then relate the force terms due to the outside singularities to the induced 
velocities inside the body and express the forces in terms of the relative 
induced velocities and strengths of singularities inside of the body (induced 
velocity method). The induced velocity method is then extended to the case 
of multibody flow. Both the singularity velocity method and the induced 
velocity method, for discrete singularities, are finally extended to the case of 
continuously distributed vortices, sources and doublets. 

In section 3, we use three problems to show how to use or to validate 
the force formulas. The first problem is for circular cylinders for which the 
method of images can be used to obtain the singularities. Notably, we will 
consider the interaction of two circular cylinders with given bound vortices, 
for which Crowdy (2006) gives exact solution. The second problem is the 
drag for Karman vortex street, a very difficulty problem since the drag is not 
related to the shape of the body. The last problem is the interaction of a 
free vortex with an airfoil, including the well known example of impulsively 
starting flow with vortex shedding. Finally we will consider an application, 
a bound vortex above the middle point of a flat plate. 

A short summary, with emphasis on the new features of the present work 
and remaining work to be done, is provided in section 4. 



2 Force formulas in various forms 

In this section, we first state the flow fleld to be considered. Then we use a 
momentum approach based on a suitably designed control volume to obtain 
a force formula. Then this force formula is rewritten in a form such that the 
forces are only related to flow properties inside the body. Finally the force 
formulas will be extended to multibody flow and to flows with continuous 
distribution of vortices, sources and doublets. Important remarks, including 
relation and difference to other theories and the treatment of vortex produc- 
tion forces, will also be provided. 

2.1 Description of the flow fleld 

Consider a body with infinite span, immersed in an incompressible two- 
dimensional flow at constant density p. The freestream velocity \4o is 
assumed horizontal. The local flow fleld is supposed to be generated by 
vortices, sources, doublets and body acceleration and rotation, in a way that 
the total velocity of the flow can be obtained by a linear superposition of 
the induced velocities due to these factors. Singularities, including point 
vortices, sources and doublets, are assumed to be either inside of the body 
(called inner ones) or outside of the body (called outer ones). Each of the 
outer singularities will be assumed to be at a flnite distance to the body. 

The sum of the strengths of the inner vortices is equal to Ft = /^^ {udx + vdy), 
which is just the circulation of the bound vortex, the closed curve dA is along 
the body with an anticlockwise path, so that a clockwise circulation has a 
negative sign. We note that even when there is vortex production, the con- 
servation of total circulation holds 

i 

A singularity, located at {xi,yi) but generally moving at the velocity 
(^, -^), will be either a point vortex of strength Fj, a source of strength 
mj, or a doublet of strength /Xj. 

The (fluid) velocity induced at {x,y) by a point vortex (i) at {xi,yi) is 

uii)(x y) = ^ = ^ = r, y-Vi 

y(^)(a;,y) = -^ = ^ = ^S^ ^^^ 

\ T !^ / dx ov 2-K r- 



where ijji = — |^ Inrj with Tj = ^/(a; — Xj)^ + (y — ViY is the stream function 
and (pi = -^ is the velocity potential. The angle Oi is defined such that 
X — Xi = Vi cos 9i and y — yi = rj sin 6*4. It should be emphasized that this 
induced velocity is independent of the velocity of the vortex. 

A doublet can be treated equivalently as a vortex pair (Thomson 1968,p361), 
see at the end of section 2.3. Thus we first assume the doublets have been 
transformed into vortices and just derive forces due to vortices and sources, 
then the explicit infiuence due to doublets will be derived directly from the 
vortex based forces (section 2.3). 

Pure source (sink) singularities are in fact not required since we only 
consider closed bodies, for which we may always use a number of source 
doublets to represent pure sources. But for completeness we will also consider 
the existence of sources. The fiow field due to a point source (z) of strength 
rrii is 



dx dy 2tt 



dy dx 2tt rj /<j\ 



Since the functional forms for ([3]) and ([2]) are similar, the forces due to sources 
can be similarly obtained as for vortices. 

Now consider body generated fiow, for a body {A) rotating at the angular 
speed n (positive if anticlockwise) around the point (xq, yo) which translates 
in addition at velocity (f/, V) in a fiow already with a free stream velocity 
Vcxi- The fiow potential and stream function due to body translation and 
rotation may be decomposed as 

(Pb = U(t)u + V(t)v + ^<Pn, i)b = Ui)u + Vil)v + ^i^n (4) 

where (pUy'Pv, 4>n and ipu, i^v, i^n are the so called normalized potentials and 
stream functions, generated by the body translating and rotating at unitary 
speed. The forces due to this will be related to added mass effects. 

2.2 Singularity velocity method 

Now we use a momentum balance approach, based on the control volumes 
(vertical control volume and horizontal control volume) defined in Fig.l, to 
study the forces. The use of two types of control volumes is just for simpli- 
fication of some algebraic operations during the derivation of forces. 

We assume the boundaries of the control volume to be far enough away 
from the body and singularities, so that the momentum balance approach 
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Figure 1: Ccontrol volumes: the solid box ABCD defines a vertical control volume 
{Xi, X2) X (—00, cxd) and the dashed box abed defines a horizontal control volume 
(-00,00) X (Fi,F2)- 



used here will be linear and therefore the contributions by various singular- 
ities (when the body is regarded fixed) and by body generated flow (when 
the body is accelerating and rotating) can be decomposed. 

A) forces due to vortices. For lift, we use the vertical control volume 
for momentum balance. The body is subjected to a lift force (L^), due to 
vortices, so that the fluid in the control volume is subjected to a force of 
equal magnitude but with an opposite direction, and this force is balanced 
by the momentum flux across the left and right boundaries {x = Xi, x = X2) 
and the time variation of the momentum inside the control volume excluding 
the body, i.e.. 



(i) 



^ J X = X\ ^ J X=X2 j "''^ 

The last two terms on the right hand side represent the momentum change 
in the control volume excluding the region occupied by the body. Hence if 



dm. 



dt 



is deflned for the whole space in the control volume, i.e. 



— pv^ 'ayax 



dt dt 



Xi—00 

then. Lav represents the momentum variation rate of the flctitious fluid inside 
the body. The speciflc role of Lav will be further discussed in the end of this 
subsection. 

To flnd the explicit form of the integrals involved in ([5]), we use the 
identity 



;dy 



TTC 



which holds for any set of parameters 6, c, d independent of y. Hence 

00 

yj ^00^ IXi y 2 \Xi-Xi\ 2f^^°°^'- 
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/ / ,.«dydx = p^ll r./ tzif + r J ^(^ 



Xi — 00 



Xi Xi 



dx 



pMt (-^ (^i - ^1) r^ + TT (X2 - Xi) Ti) 

dt 
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-pr.fe)+p^(vr(x, + xo)t 



Inserting these formulas into (^ we get 



dt 



L. = -pJ2 (^-r. - ^^^) + La. + p^ (vr(X2 + Xi)) J2- 
Using ([1]) to eliminate the last term on the right hand side we obtain 

L. = -pJ2 (^-r. - ^ffj^] + L„. (6) 

Similarly, by using the horizontal control volume, we obtain the drag force 
formula 

i 

where Day represents the change of x-momentum inside the body. 

B) forces due to sources. As remarked in the last subsection, the func- 
tional form of the velocity components u and v due to a point source is the 
same as that of v and —u for a vortex. Hence we can use the results of 
vortices and directly write down the force formulas for point sources as 

Lst = pJ2 \t + ^-' ^st = pJ2 \l + ^- (8) 

i i 

and Las and Das represent the momentum change inside the body due to 
sources. 

C) momentum change inside the body. Now consider the momentum 
changes inside the body due to vortices {Lav,Dav) and sources {Las, Das)- 
With La = Lav + Las and Da = Dav + Das we may write 

with the body fixed since the role due to accelerating translation and rotation 
will be treated separately below. 
In appendix A, we will prove that 

L, = 0,D„ = (10) 



D) body acceleration and rotation, added mass effect. The forces {Ladd, Dadd) 
due to body acceleration and rotation have been well studied in the past us- 
ing either the kinetic energy method (cf Lanib(1932)) or the unsteady Blasius 
equation (cf. Wu, Yang &Young (2012)), and have been shown, by Wu, Yang 
&Young (2012), as 

Ladd = — ^ {UAuu + VAuv + ^Auu) /-, -| \ 

Dadd = ^di (L^^uv + VByy + ^B^q) 

where Auu,Au^,Auq, and i?u„,i?„^, i?„n are added mass coefficients. The gen- 
eral method for computing added mass coefficients can be found in Lamb 
(1932). 

E) Summary. Summing the force components defined in (in]),(IZl), dHD and 
( ITOj) . we obtain the lift and drag forces as 

L = -pE (K»r, - tei) + p^^i^ + Ladd 

* i fl2l 

^ — PZ^ dt "^ PZ^ dt "^ ^^dd 



I 



Here each sum is performed over all the inner and outer singularities. The 
relation of flT2|) with other known theories will be discussed in section 2.6 
(Remark 2.1). 

2.3 Induced velocity method 

In (TT2l) . the forces are related to the speeds of the singularities. Now we 
replace the speeds of the free singularities in terms of the induced velocities. 
For convenience, we use the condition ([T]) to make the term — pX^i^Tj 
disappeared in flT2|) and rewrite f fT2|) as 



^ = pE^r,+pEt-. + ^. 



J,OU J,OU 



j,ou j,ou 

Here (L^, D^) includes all the rest terms 



(13) 



^r - P2^ [ dt + dt + P2^ \^J dt +yi dt + ^'^dd 

k,in ^ ' j,ou ^ ^ 

n — _nY^ ( djykTk) _ d(mkXk) \ _ v^ ( ^, '£± _ ^r ^tuA _l_ 71 

^r — PZ_, \ dt dt I Pf^ l^i dt -^0 dt I ^ • 



10 



The symbol ^ means summation over all the vortices and sources outside 

j,ou 

of the body, while ^ means for those inside the body. 

k,in 

The velocity (-^, -^) for any free singularities involved in (TT3|) is due to 
freestream convection and induction by all the inner and outer singularities 
except itself, i.e.. 



-^ = Voo- Yl O^kYjk - rrikXjk) - Yl i^i^ji - ^i^ji] 

k,in l,ljtj,ou 

k,in l,l^j,ou 



Here 



Y "^i "^^ -y yj y^ 



with dL = {xj — XkY + hjj ~ ykY- Inserting this into (IT5|1 to replace the 



factors dxj/dt and di/j/dt, we may write 

r _L r _L r . . 

(15) 



L — Lf) + Lif + Lff + Lt + Lp + Ladd 

where the various components on the right hand sides are given and discussed 
below. 

a) The component {Lf,,Db) is 



Lb = pV^J2rj, Db = 



J,OU 



Since ^b = — J2 ^j is equal to the circulation of the total bound vortices. 



3,ou 



we have Lf, = — pVooFft and Df, = 0. This is just the basic force given by the 
Kutta-Joukowski theorem. 

b) The force {Lij, Dif) defined as 



^i/ = -pE E (^kYjk - rUkXjk) Tj - E (^kXjk + mkYjk) rrij 

j,ou \k,in k,in 

Dif = -pYj Yj {^kXjk + mkYjk) T^P + Y. (rfc^^-fc - mkXjk) ruj 

j,ou \k,in k,in 
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is due to the interaction between the inner singularities and outer singulari- 
ties. Putting those terms with a factor F^ (and similarly m^) together and 
then exchanging the order of the double sum, as 

Lif = -P E TfcXl (Xjk^j - ^jk^rij) + pJ2^kYl i^jk^j + ^jk^rij) 

k,in j,ou k,in j,ou 

Dif = -pYj^kYj i^jk^j + ^jk^j) - pYj'f^kY. yjk^T - ^jklTT-j 
k^in j:OU k^in j:OU ^ 

we obtain 

^if = "/^ 5Z ^''^'^ ^ P'l^ '^kVk, Dif = p^ VkTk - P^ Ukirik (16) 

k,in k,in k,in k,in 

where {uk,Vk), defined as 

^k = Yl (^Jfc^i " ^Jkruj) , Vk = -Y^ i^jk^j + Yjkrrij) (17) 

j,ou j,ou 

is the fluid velocity, at the location of the inner singularity {k), induced by 
all the outside singularities. 

c) The force {Lff,Dff), defined as 

Lff = pE E O^iYji - miXji) {rrij - Tj) 

j,OU l,ly^j,OU 

Dff = -pY. E (r;X,/ + m^F,OK + r,) 

j,OU l,ly^j,OU 

is due to the mutual interaction between the free singularities. It is obvious 
that the contributions to this force by each pair of j, I with j ^ I mutually 
cancel and thus 

^ff = ^ff = 
Hence the force due to mutual interaction between the free singularities does 
not contribute to forces. 

d) The force component {Lt,Dt) defined as 

r _ sr^fd{xiTi) d{ykmk)\ _ ^^fd{ykTk) d{xkmk)\ 

^'-p2^ \-iit d^j ^^'--p2^ \r^r~ + ~^^) 

k,in ^ ' k,in 

is due to the motion and production of strengths of the inner singularities. 

12 



e) The force component {Lp, Dp) defined as 

j,ou ^ ' j,ou ^ ' 

is due to production of vortices and sources outside of the body. If there are 

no vortex production, then —^ = for each free vortex and Lp = Dp = 0. 

f ) In the above derivation, the doublets have been grouped into vortices, 
since each doublet can be represented by a vortex pair. Now we make the 
force contribution due to doublets (inside the body) in an explicit form. As 
shown in section 2.1, each doublet of strength /ij and at position (xj,?/,) can 
be considered as a vortex pair of strength ^Ti at (xj,?/j ± e) with e — )■ 
and 2eri = fii. Apply flT6|) to the corresponding vortex pairs yields, for 

e -^ 0, a force component {L^, Du) with L^ = P Z]!=i ^^2£:rj^\ D^ = 

-^-^dui sr-^ dvi 

i,in i,in 

The various force components above will be put in compact form in section 
2.4. Before doing this we would like to remark that the above analysis appears 
to have given a way to interpret the physical origin of each force component. 
This is not seen elsewhere according to the knowledge of the present authors. 
The force due to vortex production will be further discussed in sections 2.4- 
2.6. 

2.4 Summary of the induced velocity method and multi- 
body extension 

Now assume there are vortices (not including the doublets now), sources and 
doublets inside the body, and outside of the body there are a number of free 
vortices and sources. Inserting the force components defined in items a)-e) 
in section 2.3 into (II 5p we obtain the force formulas below 

L = —pVoo^b + Lind + Lt + Lp + Ladd /-, q\ 

D = Dind + Dt + Dp + Dadd ^ ' 
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Here {Und^Dind), defined as 



i,in i,in i,in {nr\\ 

Dind = E P'^i^i - E P'^i'^i - PJ2 % Pi 

i,in i,in i,in 

is due to the induced velocity effect at the inner singularities, and (Mj,fj), 
defined by (fT7|) and rewritten here as 

_ , , _^ji.yi-yj) , ^A^i-^i)\ „_\^ /^rj(xi-Xj) ^ mj{yi-yj) 

Uj, 






27r4 27r4 7 ' ' ^ V Stt^^^ 27r4 

(21) 
denotes the fluid velocities induced at the location of the inner singularities 
by all the outside singularities, not including induction by other inner sin- 
gularities. The force component {Lt,Dt), due to motion and production of 
inner singularities, is defined by 

i,in ^ ' i,m ^ 

(22) 
Finally the component {Lp,Dp), defined by ( ITSj) . is due to production of 
singularities outside of the body. Due to Kelvin theorem of conservation of 
circulation, we have -^ = once the vortex is produced and is moving freely. 
Hence it remains only those just in production. Generally, vortices will be 
produced at some geometric singularities, such as the trailing edge of a body. 
Moreover, we do not consider the possible case of source production outside 
of the body. Then 

Lp = PZ^^;, — , Dp = -P2_,ys-^ (23) 

s s 

Here the summation is performed over the points on the surface of the body 
where we have vortex production at rate dTg/dt. The usefulness of this 
term will become clear in for instance the treatment of Wagner problem with 
vortex production. 

The similarity and essential difference of flT9|) comparing to the force 
formula of Wu, Yang &Young (2012) will be discussed in section 2.6 (Remark 
2.3). 

14 



Now we discuss the extension of the force formula fll9p to the case of 
multiple bodies. This force formula has been obtained for a single body 
without using pressure integration. If the pressure p is used to integrate the 
force, a.s L = § pdx, D = — § pdy, we of course should have the same forces 

dA dA 

as given by ( IT9|) . Now remark that the forces in the form of (fT9|) only depend 
on the induced velocities inside the body and the motion and production 
of singularities inside and on the body, and that the flow pattern (such as 
induced velocity and their derivatives) inside the body is the same whether 
the outside singularities are free ones or bound ones (as created by another 
body). We thus have the same force formula if there are outside bodies, 
provided the flow induced by the outside bodies be represented by a flow 
induced by equivalent singularities. 

Consider, for the case of multiple bodies (namely body A, B, C, ■ ■ ■), 
the force formulas for the body A with contours dA. Then, according to the 
above remark, the force formula for body A is 







(24) 



where FJ, is the circulation around body A, 



i,A i,A i,A 

^Ind = E P'"^'^i - E P^i^t + P E ^Pi 

i,A i,A i,A 

is the induced velocity effect, with summation performed over all singularities 
inside body A, {ui,Vi) is the velocity at {xi,yi) induced by all the outside 
singularities and bodies. The unsteady term {L\ \Di ), now defined by 

r(^)_ >v^ ( d{xiTi) d{yimi) \ ^(a) _ s^ ( d{yiTi) ^ djximi 

^* -pl^y-^t d^ ' ' --p^y-^^ + ^dT 

i,A ^ ' i,A ^ 

is due to the motion and production of singularities inside body A. Finally, 

s,A s,A 

is due to the vortex production (at for instance geometric singularities) on 
the surface of body A. 

The force formulas for bodies -B,C, ■ • • can be similarly defined. 
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2.5 Force formula for distributed sources 

Let oj = ijj{x,y,t), m = m{x,y,t) and fi = ii{x,y) be the continuous distri- 
butions of vortices, sources and doublets, and let {u,v) = (u{x,y),v{x,y)) be 
the velocity inside the body and induced by all the outside vortices, sources 
and doublets, free or body generated. With the help of Dirac delta function 
we may transform the force formulasabove into integral forms. 

A) Singularity velocity method. For the force formula (fT2l) . the integral 
form is 

L = PTtI Id^ ujxdxdy + pfj J^^ mydxdy + Km ,^^. 

D = -pit I Id^ ujydxdy + pfj J^^ mxdxdy + Dadd 

where D^ is the region occupied by both the fluid and solid. Here we have 
used / /^ udxdy = 0. 

B) induced velocity method. Consider only multibody problems since the 
single body problem is only one such a special case. For body A, the integral 
form for (^ and ([24D is 

La = -pVoo J Xi ujdxdy - p J Jj^ (uu + mv - |^/ij dxdy + lJ^^ + lJ,^^ + l[^1 

Da = pJL (^^ -^^ + i-") dxdy + dJ^^ + 4^' + D[f, 

(26) 
where 

(27) 



and 



4-''=''E-.§.'Dr' = -pE!"§ (28) 



dt' ' • '^" dt 

s.A s,A 



Both 0251) and (I26p will be validated against the problem of an impulsively 
starting plate (section 3.3), where we have vortex production. We should 
emphasize that the influence due to vortex production inside the body is 
embedded in (12] 



2.6 General remarks 

Now we provide some general remarks about the force formulas derived above, 
including the connections to known theories and the new features. 
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Remark 2.1. The formula f ll2p (singularity velocity method) is rather 
general, including the contribution from both the inner and outer singulari- 
ties. The way to express the forces and the way to obtain this formula are 
new. Its integral form (1251) is the same as the integral approach of Wu (1981) 
for inviscid two dimensional flow, except we have an additional term due to 
sources. When Fj is constant and when the body is fixed, the force formula 
flT^ simplifies as 

i-E^(^-§)r..^--E^|r. (29) 

In the special case that the outside vortices, for instance the starting vortices, 
move at the freestream speed and the internal vortices are fixed, as in the 
case of steady flow, we recover the classical KJ theorem from fl29l) 

L = -pVooTb, D = (KJ theory) 

Here Tb = X^Ti- In order to recover the Kutta Joukowski theorem, the far 

in 

away starting vortex must be taken into account by the integral approaches 
of Wu (1981,p438) and Howe (1995,p416), even for steady flow. When the 
reference frame is such that the far flow field is still, an artificially built 
image distribution inside the body is needed to recover the Kutta Joukowski 
theorem in the approach of Saffman (1992, p48). The forces expressed in 
the form of fl29|) or flT2|) mean that it is instead the relative velocity of each 
vortex that determines the force. For instance, for the free vortex i, the force 
components are proportional to its relative velocity (\4o — ^; ~l^)- Hence 
the way to express the forces here is frame independent. 

Remark 2.2. The forces due to vortex formation, for instance in the form 
of (l23ll . is strange due to the dependence on the position Xg and z/g. This 
would mean that the magnitude of the forces depend on the choice of the 
reference frame. This is in fact not so since in real problems the vortices 
always produce in pair due to conservation of vorticity. This means that the 
production of one vortex of circulation F„ at Xn (which may be some point 
near the trailing edge of an airfoil) is at the consequence of the production 
of another vortex of circulation — F„ at a point x„ (which may be a point 
inside the body close to the trailing edge) close to x„. Since the momentum 
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due to this vortex pair is 

X2 00 

= ^p(y(Xi+X2-2a;„)-^(Xi + X2-2J„)') 

A: ^ '' 

n 

the forces due to this production, when motion is excluded, are thus 

Lp = -^p^{Xn-Xn),Dp = ^p-^{yn-yn) (30) 

n n 

Hence the magnitude of the force due to vortex production is frame inde- 
pendent. The way to treat the vortex production outside the body in the 
way of expression ( l28l) and that inside the body embedded in ( 127|) is very 
convenient, see section 3.3 for the Wagner problem with vortex production. 

Remark 2.3 Compared to Wu, Yang &Young (2012), where the force 
formula has been obtained directly through the unsteady Blasius equation 
(though suitable only for irrotational flow so that it does not apply to the case 
when vortices are produced on the surface of the body), the force formulas 
f lT9|) and fl24|) based on the induced velocity method are more general since 
here we include the role of bound vortices and vortex production. Moreover, 
the induced velocity in the force formula of Wu, Yang &Young (2012) is 
due to all the singularities (including the inner ones) while the present one 
this induced velocity is only due to outside singularities (and bodies). As 
remarked by them, the contributions to induced velocity effect from any 
pair of interior singularities cancel out and only the contributions from the 
external vortices remain. Thus both approaches yield the same force for 
induced velocity effect. 

Remark 2.4 The force formula (1261) is in fact some new form of integral 
approaches, valid for multiple bodies. The past force formulas based on inte- 
gral approaches involve volume (and sometimes boundary) integrals defined 
either in the entire space or in the fiuid regime, and requiring the use of auxil- 
iary potential functions for multibody force decomposition. The present one 
uses integral only defined inside the actual body and has the advantage that 



it does need auxiliary potential functions. As shown in section 3.3, it will be 
very convenient to be used for airfoil problems where the internal singularity 
distribution can be found through standard methods. 

3 Validation and application 

The force formulas given in section 2 rely on the knowledge of the circulation 
of the bound vortex and its time variation, position and speed of the vortices, 
sources and doublets, in discrete form or distributed form. In this section we 
give several examples to demonstrate the application of the force formulas. 
In section 3.1 we apply the induced velocity method to circular cylinders 
for which the singularities can be determined by the method of images. In 
section 3.2 we use the singularity velocity method to study the force for 
Karman vortex street, which involves an infinite number of discrete vortices 
and the force formula need be specially adapted to this case. In section 
3.3 we use the integral form of both the singularity velocity method and the 
induced velocity method to study the lift of a thin airfoil with vortex shedding 
or interacting with another airfoil presently represented by a lumped vortex. 

3.1 Problems of circular cylinder 

Here we first simplify the induced velocity method for the case of circular 
cylinders, then we apply the results to study one cylinder with a pair of out- 
side standing vortices and a source doublet, and the problem of two circular 
cylinders with given circulation. All the problems have known solutions so 
that they are used to validate the present force formulas. 

3.1.1 Simplified force formula for a circular cylinder 

For each vortex j of circulation T^ outside of a circular cylinder, there is one 
image vortex of circulation T ° = T- at the origin, and one of circulation 
Tj = —Tj at the inverse point {Xj ,yj ). An outside source doublet at 
{Xj ,yj ) and with strength o"j, has an image doublet of strength 

(m) O, 

^3 = J2^i 
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at its inverse point, where fj is the distance of the outside doublet to the 
body center. The force formula (I19p applied here yields the following force 
decomposition 

L = Lb + Lh + Lo + Lm + L^ + Lf + Ln 



with 



Lb 
Lb-- 



-pVocTb, Db 







-p Yl UiT] ', Db = p^ViT 



(b) 



^(o) 



Lo = -pY1 ^i'^i , Do = pYI ^i^ 



(o) 



Lr 



L, 



L, 



-p J2 uiT 



(m) 



; -'-^m 



pE^«r 



(m) 



(31) 



i,in 

Edui a? a 
— 
j,ou 

0, Df 



dy 



W 





D. 



EdvH 
— 



j,ou 



dy /? 



Lf = u, lyf 

Here {Lb,Db) is the basic bound vortex force, the components {Lb,Db), 
{Lo, Do) and {Lm, Dm) are due to induced velocities at the locations of bound 
vortices, image vortices at body center and image vortices at inverse point, 
respectively. The component (L^, D^) is due to the induced velocity gradient 
at the location of inner real doublets, and the component {L„, D^) is due to 
the images (at the inverse point) of the outside doublets. The velocity (uj, Vi) 
is the induced fluid velocity (induced by all the outside vortices and source 
doublets) relative to the speed of the internal singularity [i). Finally the 
component (L„, D^) with Ln = -Y.nP 



dt 



Xri. •^i 



is due to vortex production, (x„,|/„) and (x„,y,i 
vortex pair with circulation production. 



denote the position of the 



3.1.2 Standing vortex pair behind a circular cylinder 

It is well known that at moderate Reynolds numbers, the flow around a 
circular cylinder involves two standing, oppositely rotating vortices behind its 
cylinder. An inviscid model for this consists of two equal and opposite point 
vortices, of circulation F > and — F < 0, standing symmetrically behind 
the cylinder (Saffman 1992, p42, Milne-Thomson 1968,p370), at the positions 
x^ = {rj cos9f,rfsm6f) and lr_ = {rj cos9f,—rfsm6f), respectively. 
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The image vortices at the inverse points are respectively at Xj^ = 
{j-cos9f,^sm9f) and 1^_ = {^ cos9f,—^sm9f). On the Foppl hne 
(see for instance Saffman 1992,p.43) defined by 

(rj -a^Y = 4rj sin^ {9f) or 4r) - (rj - a')' = 4rj cos^ {9f) , 

the two vortices, though under the convection by stream fiow and under 
induction by the vortices (including images) and source doublet of strength 
H = 27ra^Voo, remain stationary if the circulation is given by 

r± = ^r, r = AnV^rj sin {9f) | 1 - ^ 

It is well known that for this case the drag vanishes, either by direct calcula- 
tion of pressure on the cylinder or by considering the vortex pair as a source 
doublet at far enough distance (Saffman 1992, p43). Now we will check if we 
recover this conclusion by the force formula in terms of the induced velocity. 
The induced velocities v^ and v_ at the two inverse points and the in- 
duced velocity v^°' and its derivative ^^ at the center of the cylinder are 
found to be 

(m) _ a^VooSm(2ej) („) _ a^Kx, sin(2e^) 

V+ — Z^ ,V_ — -2 

' dy vr r^ 

The drag force due to the induced velocities at the two image points is, 
due to the force relation in (13T|) . 

D^ = pT^r^v^r^ + pT^-V_-^ = 2pra Voo^sin {29 f) 

and the drag due to the doublet at the body center is, according to the fifth 
relation in fl3T]) . 



D = -pji 



at;(°)(0,0) 2a2p14orsin2^ 



/ 



M ^^^^ Q 2 



dy rj 

Hence Dm cancels D^, and the total drag obtained by the present method 
vanishes. 
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3.1.3 A doublet outside of a circular cylinder 

Consider a doublet of strength a at y = f, then the velocity induced by this 
doublet on the line x = is ^'•'^•'(0, y) = 2^(1^- )2 ^° ^'^^^ ^^^ derivatives at the 
center and inverse point are 

-(0,0) = — -,^— (0,y„ 



dy ' ' ' nP' dy ' '^'"^ nif-a^fY 

According to the fifth and sixth relations in flHT]) . if there is a doublet of 
strength fi at the center of the cylinder, there is a lift force given by 

and there is an outside doublet, there is a lift due to the image of this doublet 

The latter is the same as that given by the Blasius theorem or by Lagally 
theorem (cf. Milne-Thomson, 1992,p232). This is a force which points to the 
doublet. 

3.1.4 Two circular cylinders with circulation 

Crowdy (2006) gives a general theory which permits the calculation of lift 
for a finite number of staggered cylinders with bound vortex. Here we apply 
fl3T|) to obtain the lift for his example of two vertically aligned cylinders 
(see Fig. 2), both of radius a = 1/2 and of given circulation Fi = r2 = F 
(F = —5,-2,0), immersed in a uniform stream \4o = 1 and placed at a 
distance h > 1 between their centers. The bound vortex associated with 
one cylinder has an image pair of two counter rotating vortices in the other, 
at the center and inverse point respectively. An image pair in one cylinder 
also has two counter rotating image vortices in the corresponding inverse 
points of the other cylinder, plus two cancelling images at the center of the 
other cylinder cancel. Thus in each cylinder, there are an infinite number of 
such inverse points, with counter rotating vortices of equal strength between 
any two adjacent inverse points, with distance becoming closer for newer 
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V. 





V, 





Figure 2: Staggered two cylinder problem of Crowdy (2006). Left: staggered 
cylinders with bound vortices. Right: vortex system with image vortices for 
the equivalent inverse point model. 

generated images. If, for the ith cylinder, the distance of the kth inverse 
point to the center of this cylinder is denoted as /i^ , then 



h^^h 



and hu I h 



h-h 



(0 

k-l, 



a for k > 1 



(32) 



For an exact solution with (EI]), we have to work with such infinite 
number of inverse points. Here we instead use an approximate method based 
on the remark that the inverse points in each cylinder are distributed in a 
narrow region. To see this, let hj^ = g^^^ for k — > oo, then by ( !32|) we have 



*) h- 



h-g 



{«) 



which can be solved to give 
g^^ =-{h- Vh^ - iaA 



(33) 



and one can verify that a'^/h < h^ 



''k ^ 9^^^ f°^ k > 1. It can be further 
verified that the region between y = a'^/h and y 
for h close to 1. 



g^"^^ is very narrow even 
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Hence within the framework of approximate solution we may merge all 
the inverse points into an equivalent one, with a distance he to the center 
of the corresponding cylinder satisfying a^/h < he < g^^^ ■ Moreover, since 
these inverse points are denser close to y = (?'■*•' , we just set he = a^/h. 
With a = |, we have he = he = 0.5'^/h. Let hee be the distance between 
the equivalent inverse points of the two cylinders, and heo be the distance 
between the equivalent inverse point of one cylinder to the center of the other. 
Then 

hee = h- 2he, heo = h - he 

With the approximation of equivalent inverse point, and considering the 
conservation of circulation, the vortex system in each cylinder can be simpli- 
fied in the following way. 

For the 1st cylinder, there is one given bound vortex of circulation T[ = 
Fi and one image vortex of circulation T° = F2 at the center of this cylinder, 
plus one equivalent image vortex of circulation F^™ = — F2 at the equivalent 

(i) 

inverse point y = hi . For the 2nd cylinder, this can be similarly defined. 
Moreover, each cylinder has a doublet of strength /ii = /X2 = /^ = 27ra^Voo = 
|. There are also an infinite number of image doublets at the inverse points. 
When the equivalent inverse point is applied to the doublet, the strength of 
the doublet at this inverse point is ^i = ^2 = Z^*"'^'' ~ f^'^'^/^lo- 

The velocities induced at the center and at the equivalent inverse points 
of the 1st and 2nd cylinders by the bound and image vortices and doublet of 
the 2nd and 1st cylinders are thus 

,(°) _ „,(o) _ 2r I -r I M I m'"' 



,,(™) _ _,,(™) _ 2r , -r , fi , j/=) 

"1 — "9 — o„!, "T n^h ' n„l,2 ' 



'1 "2 2nheo ~^ 2TThee ~^ ^Trhl^ "^ 2-nhl^ 

The derivatives of the corresponding induced velocities are 

9n^ ^ _a4^ ^ ^r^ -r , ^ , n(<=) 

dy dy 2-Kh'^ """ 2-Khl^ "*" irh^ "*" -nhl^ 

dy dy 2-Khl, "^ 2-Khl, ^ nh^, ^ tt/iL 

According to (!3T]) . the force formula for the lower airfoil is L^^^ = L^ +Ll + 
L^}^ + L« + L« with L« = -pl/ooFi, Li^) = -p«S°)Fi, L« = -pu'^^T^\ 
r(l) Mt-M r(l) ^^i"' r(l) ^"i"' (e) 
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With p = \4o = 1, Ti = Fg = — r2 = r and yUi = /i = ^vr, we have 

ay aj/ 

where the induced velocities u" ,Ui and their derivatives du[ /dy, du{ /dy 
have been given above. Similarly for the upper airfoil, 

dy dy 

The results, compared to Crowdy (2006), are displayed in Fig 3. We 
remark that the agreement is acceptable even when h is short and despite 
the use of equivalent inverse point to merge all the inverse points. The short 
distance behavior, that is, there is an attraction force when the cylinders are 
close, has been discussed by Crowdy. Here it is found that this is due to the 
influence of the real and image doublets and the induced velocity gradient, 
described by L^ and L„ in (13T]) . 

3.2 Karman vortex street 

The problem of the Karman vortex street behind a bluff body (Fig.4) is 
rather special since the shape of the body is unknown. Hence we can only 
use the singularity velocity method to study the drag. The Karman vortex 
street is a double row of staggered and counter-rotating vortices of strength 
F > and F < and moving horizontally at speed 

Fvr Tib 

V = — tanh — 
a a 

Here b is the vertical separation distance between these two rows and a is 

the horizontal separation distance between adjacent two vortices in each row 

(see for instance Milne-Thomson 1968, p377). The period of vortex shedding 

in each row is thus 

a 

According to (130!) . the drag, averaged over r, due to the shedding of new 
vortex pair, separated at a distance yn — yn = b, is 
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Figure 3: Comparison between the present result (solid lines) with Crowdy 
(2006, Fig.4) (dashed lines), for F = (lower), F = —2 (middle) and F = —5 
(upper). Dashed hues: Crowdy. 



A, 



P- 



1 f d Vn iVn -yn)dt _ p_ 

dt 



d{Tn{yn-yn)) 



pbT 



r Jo «r -r JO 

When the expressions for V and r are used, we obtain 

pbT {V^ -V) _ pbT ( _ Ftt 

a 



A, 



Voo tanh — 



a 



This is the well-known formula for the unsteady part of the drag, which has 
been otherwise obtained by sophisticate Blasius approach (Milne-Thomosn 
1968, p382), impulse approach (SafFman 1992,p 137) or integral approach 
(Howe, 1995, p417). The way to obtain this force using the present approach 
appears to be much easier. 

Remark also that, apart from the drag component Du due to shedding 
of new vortices, there is also another component Ds due to quasi steady 
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Figure 4: Control volume for the Karman Vortex street problem. A outline 
X = Xd separates the flow regime into an upstream part and a downstream 
part (D). 



flow formed by two rows of periodically counter-rotating vortices spreading 
infinitely into the downstream direction. No existing force formula, including 
the present formula f lT2|) . can be directly applied to this case. One then needs 
to construct a special downstream boundary and relates this component of 
drag to the momentum flux across this boundary. For details, see Milne- 
Thomosn(1968,p382), (Saffman 1992,p 137) and Howe (1995,p417), where 
different ways for setting the downstream boundary are used. Let x = Xd 
be such a boundary downstream of the body and intersecting the vortex 
street in its uniform region. Let {uk, Vk) be the induced velocity by the 
vortex street. Then, in Appendix B, we shall prove that, the use of (1T2|) 
leads to 



D.. 



2P 



Ui 



I) dy 



(34) 



which is exactly the same as obtained by Howe (1995, Eq.(5.11)) through an 
integral approach. Howe then shows that fl5^ yields 



D. 



2a? 



tanh — 

TT a 



The total drag D = 
D 



Du + Dg is thus 
phV {V^ - V) 



+ 



2a2 



a Tcb 
tanh — 

vr a 
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3.3 Interaction of a thin airfoil with outside vortices 

Consider a thin airfoil without thickness, for which the camber hne is defined 

by 

y = ya{x) = -ax + f{x), < x < ca (35) 

Here a is the angle of attack and f{x) is small with respect to ca- The 
upstream inflow is assumed horizontal. First we need to find the distribution 
of vortices. Then we apply both fl^ (singularity velocity method) and (12B]) 
(induced velocity method) to obtain the force formulas. The results will be 
validated through the problem of Wagner with vortex shedding. Finally the 
interaction of the airfoil with another one, represented by a bound vortex, is 
studied. 

3.3.1 Method of solution 

To solve this problem we assume a distribution of vorticity 7(x) on the airfoil 
so that the velocity at any point on the airfoil is 



J 2tt {{x - + [y ~ yaiOr) 



CA CA 

J 2n {{x - 0^ + [y ~ yaiOr) J 27r(x - 



where {uv{x),Vv{x)) is the velocity on the airfoil induced by vortices outside 
of the airfoil and is the only correction to the standard method for thin 
airfoil theory (Anderson 2010). Substituting the expressions of ua{x) and 
Va{x) into the boundary condition 

va , df{x) 

— —a + 



ua dx 

we obtain the equation for 7(x), 

CA 





where the function 



[ iJQd^ ,r (df\x) \ 



Uv[x) fdf{x) \ Vv{x) 



g[x) = ——- — a 

Vr^, \ dx 1 Kx 
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is due to the induction by the outside vortices. The standard thin airfoil 
theory is recovered if g{x) = 0. Let x = ^ {1 — cos 6*) , ^ = ^ (1 — cos (3) for 
< 9, (3 < 71, then 7(0;) can be expressed as 

7(^) = 2^00 Ao ^."'' + 5^ ^n sin (ne) (36) 



sin^ 



with 



An,n>i = -^ Jq ^ COS (n/3) (i/3 + a^ 



(37) 



where the coefficients 

an,n>l = ;J^ /J" (^^i;(0 " "^(0 (^ ~ ")) '^^^ ^^^^ ^^ 

are due to the outside vortices, and the remaining parts in Aq, An are given 
by the classical thin airfoil theory. 

Remark that in the case of a thin airfoil with thickness, we may add, in 
addition to the vorticity distribution 7(x), a distribution of source doublet 
fi{x) on the camber line. Then the boundary conditions should be defined 
for both the upper and lower surfaces of the airfoil which will provide us the 
necessary conditions to determine both 7(x) and /i(x). 

To find {ujj{x),Vv{x)) the position of each outside vortex (i) must be 
determined, through solving the equations 



/■ 



CA 

'' ^ J 2n ((x. - 0^ + {y. - yaiOf) fji, 27r ((x. - x.f + (y. - y,f) 

CA 
dy^_ f l{.i){.X^-i)di , -^^^ Vj{i){x^-Xj)di 



+ E 



* ; 2i ({x - 0- + {y- ya(0)'^) ■ f^. 2i ({x. - x,)^ + (y, - yjf) 



(39) 



3.3.2 Lift force by the singularity velocity method 

The lift force in the form (1251) can be rewritten here as 



CA ^ j-CA 

-pVoo I idx + p— / '^xdx + Lout 
/o d'l^ Jo 
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where Lout is due to the outside vortices. For point vortices of circulation Fj 






and for a distribution of outside vortices of strength k 

Lout = -pVoo / / kdxdy + P-r: kxdxdy (40) 

With 7(0;) in the form fl36l) . we get for the circulation of bound vortex 
Fb = / -fdx = ttcaVoo Mo + -Ai j (41) 

and the moment of inner vortices 



r ^xdx = ^^ r I Ao^^^ + Y An sin {nO) | (1 - cos 0) sin OdO 
Jo 2 Jo V ^i^^ ^ / 



4 V 2 

Thus the formula for the lift force is 

L = -pVooTb + Lint + Lout (42) 

where 

When the outside vortices are given and when /(^) and a are known, the 
formulas (137|) and (!38|) with ^ = ^ (1 — cos/3) yield the required coefficients 
Aq, Ai and A2. 

3.3.3 Lift force by the induced velocity method 

If we use (12^ then 

L = — pVoo / "ydx — p u'jdx + p— / ^yxdx + Lp 
Jo io dt Jo 
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where Jq'^ •yxd. 



X 



1tVooc\ 



fQ ,^«,^ — ^ (Aq + Ai — 1^2) as above, Lp = pY^ ^s^ is due 
to vortex production outside of the airfoil (that due to vortex production 
inside the body is included in the third term on the right hand side). 

Assume that for the present problem, vortex is shed only at the trailing 
edge, so that Xs = ca,^ = j^J /„„ kdxdy = - ^, and 



Lp = -pCA 



dt 



Moreover, since the vortex sheet, inducing m, is on the horizontal line down- 
stream of the plate, JJ^^ U'~^dx ~ 0, and thus 



L = -pV^T, -pi u^dx + '^ ^ ^ - [Ao + A,-^]- pcA-j^ (44) 



dt 



dt 



V« Y 
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Figure 5: Impulsively started plate with vortex shedding 



3.3.4 Wagner problem for impulsively starting plate 

Consider a flat plate of length ca at small incidence a, impulsively set into 
motion with constant velocity Vqo (Fig.5). At a distance s from the trailing 
edge, a vortex sheet of strength k{s,t) spreading over < s < Vnr,t and 

Dk 
Dt 

induced motion. 
Tb + Jo" k{s,t)ds 
k{s,t) 

k{s, t)ds + ttcaKo (Aq + -Ai] =0 (45) 



satisfying -^ = Ij + Vooff = under the assumption of negligible self- 



Moreover, since the total circulation is conserved, that is 
0, we may use pT]) to write the relation for determining 
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In Appendix C we shall prove that 

Ao + \a, = -a - —^ / °°7l - }l^A+T\ k{s, t)ds (46) 

Inserting (H6|) into (145|) . we obtain the required equation for k{s,t) 

/ ^ — A;(s,t)(is = vrcAVooa (47) 

which is exactly the same as given by the method of conformal mapping, see 
for instance Saffman (1992,plll,eq (8)). 
In Appendix C we also show that 

1 1 ^ Voo t 



Ao + Ai- -A2 = -a- —— / H{s)k{s, t)ds 
2 ttI/oo Jo 

where H{s) = J- (s — a/sca + s^) ~~ '?/^'' + ^- Using ( H7|) we may further 
write 



1 4 r°°* 4 r°°* / \ 

y4o+v4i--y42 = a — / k{s,t)ds — — 5-/ is - \/sca + s"^) k{s,t)ds 

2 ttVooCa Jo -nVooCA Jo ^ ' 

(48) 



For the singularity velocity method, the component Lint defined by (143!) 
now becomes 



A rVoat J /-Voot 

Lint = "P^A-T. / k{s,t)ds - P-f. j (■5 - a/sca + s^ j k{s,t)ds 



when (HSl) is used, and Lo„t by fllU]) becomes 



Voot J ^Voot 



Lout = -pVool k{s,t)ds + p—l k{s,t){cA + s)ds 



"'' ^0 



Inserting these into the force formula (B2|) and making use of Tb+j^ °° ^(s, t)(is 
0, we obtain 



J rVaot 

L = p— ^Jsca + s'^k{s,t)ds (49) 

dt Jq 
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For the induced velocity method, inserting ( I48p into the force formula 
[]) and remarking that J^^ ujdx ~ since the vortex sheet, inducing u, is 
on the horizontal line downstream of the plate, we get 

pcA-iT / k{s,t)ds+ p— / (s - ^/sca + s2 j k{s, t)ds j -pcA—r^ 



-v^t 



which, when Tb = — L°° k{s,t)ds is used, yields 



^0 

"Voot 



L = -pVaJ^b- P-Tf \ fs - a/sca + sA k{s,t)ds (50) 

The force formula (!50l) based on the singularity velocity method and that 
(149|) based on the induced velocity are in fact identical since it can be shown 
that 



VooTb+ jf°° skis,t)ds = (51) 

The identity (ISTjl will be proved below just for small time. 

Now we use ( l50l) to study the lift for two extreme case, t — j-oo and t — )■ 0. 
For t — J-oo, Saffman (1992,pll4) shows that 



d_ '•^-* 
di 



f s — a/sca + s^ j k{s, t)ds = 

Hence, by ([50]), 

L = L^ = -pVooTb{oo) 

where r;,(oo) = tiVooCao: is the steady state circulation. 

For small time t — )■ 0, we may write (s — a/scXT^)— J-a/sca, thus by 

(EQD, 

LK.-pVooTb+ p-T. j y/scXk{s,t)ds 

Furthermore, it can be straightforwardly verified that the solution for (H7|) 
satisfying ^ = in addition is 

k{s,tj ^ ^ ,t -)■ 

y/Voot - S 
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and thus 



Tb = - k{s,t)ds = -2Vooay/cX^/Voot, / sk{s,t)ds =VooOi^fcA 



^{Vo^tY 



/o Jo 

Due to the above expressions, flSTj) obviously holds. 
Now for small time, it holds that , F^ — )■ 0, thus 

L^p— y/scXk{s,t)ds ^ pcAVooa— ds 

dt Jo dt Jo ^/Voot - s 

Let S = y/s, then 

f , ^ ds = 2 f , ^ dS = -S^/V^t - S^+V^t arcsin -^=+C 

Hence r„^°°* ,,/! ds = ^ and therefore 

L = -TcpV^CAa 

Thus we recover the well-known result that for a flat plate impulsively 
generated, the initial lift on the plate is only one-half of the final lift Loo = 
TipV^CAd. Note also that the initial lift for airfoils with thickness has also 
been studied, see Chow& Huang (1982) and Graham (1983). 

3.3.5 A bound vortex above a fiat plate 

The lift for an airfoil with a free line vortex over the airfoil has been stud- 
ied, for instance by Saffman&Shefield (1977) using the method of conformal 
mapping. There is a condition for the vortex strength such that the free 
vortex is standing. First we apply the present theory to this case to provide 
another validation against known results. Then we replace this free vortex 
by a bound vortex, to find the force due to interaction between a flat plate 
with another airfoil represented by a lumped vortex. 

Consider an outside vortex of strength F„ > at a distance h above the 
midpoint {xm = ^) of a fiat plate of length ca = 2a parallel to a stream of 
velocity \4o (Fig-6) . The formulas (157|) and (155]) now reduce to 

I r 2 r 

Ao = -— / v^{^)dP, v4„,n>i = — r^ / v^{^)cos{nP)dP 

vrl/oo Jo ^'^oo Jo 
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t: 



Ca 



Denote h = -y—. With ft,(^) given by 

2'^ A 

r j.\ J- 1) S •^rn ^ V COS p 

27^ {^ - XmY + h^ 27rx„ cos2 /3 + 7l^ ' 

we get Aq = ^2 = and 



A, = -^ 1 - -^= = -^T{h) (^|^' + h) (52) 



where 



T{h)= f Vl + ^' + ^J (1 + ' ^^^^ 

Hence, with (136|) and (HTl) . we obtain 7(0) = 2Voo^i sin 6* and 



-2 



r^ = -T(/i) ("v/i + h' + h) r, (54) 



This condition holds independent of whether the vortex is free or bounded. 
A) Free vortex case. Apply fl5U]) to the present case, and with 7(/3) = 
2V00A1 sin/3, we obtain the velocity of the vortex (if it is free) 

^^JL -\/ _L 2T(7i)r^ ^ _ n 

''^00 ~r ; ,. U 



dt 71 C A dt 



Set ^ = 0, we obtain 



r. = -^^^ (55) 

2T{h) ^ ' 

The method of singularity approach (142|) gives a lift L = — pV^Ff, — 
P (Kxd - ^) r^, and since ^ = 0, we have 
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L = pVooir{h)l\/l + h^ + h]-l]T, (56) 

Thus if r^, is given by (!55l) and r^ is given by (15^ . then the outside vortex 
will be stationary and the lift is given by (!56|) . This result is exactly the 
same as given by the method of conformal mapping (cf. Saffmanl992,pl22). 

B) Bound vortex case. This is a multibody case, though the body above 
the plate is represented by a bound vortex. The force formula (jH]) (induced 
velocity method) applied here gives 

L = -pVooTb -pi ujdx (57) 

Jo 



Here 



U = M^,(0 



r„ h r„ h 



V 



27r (^ - x^f + h^ 2nxm cos^ P + h^ 

is the velocity in the airfoil induced by the outside bound vortex. With 
7(/3) = 2VooAi sin/3 and with Ai given by (E2]), we have 

" u,^ = Izl^il^ r ±^l^ip = -!I(5)II (58) 

TT Jo cos2/3 + /i CA 

since 

'^ —2 

sin /3 sin fidp h + 1 

Inserting ( l58l) and ( l54l) into (!57l) we obtain finally 



L = pV^Tih) [^^l + h' + h^V, + P^^^^ (59) 

Remark 3.1. The formula fl59|) with T(/i) defined by fl53l) gives the lift 
of a horizontal flat plate with another airfoil of any given circulation F^ and 
at a distance h above the middle point of the plate, provided the additional 
airfoil is simplified as a lumped vortex. 

Remark 3.2. When, in addition, the circulation T^ satisfies (155|) . then it 
is clear that ( 159|) gives exactly the same force as (156|1 . This is because that 
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when (1551) is satisfied, the vortex is standing for the free vortex case so that 
we shall have the same force as for a bound vortex. 

Define the lift coefficient as q = 2L / pV^CA and the normalized vortex 
strength as V^ = F^/V^c^, we obtain from fl59|) the lift coefficient for the fiat 
plate 



Cl 



h^ + h + 2T^] T{h)T^ 



The lift coefficient as a function of h ioi T^ = —1 and F^ = 1 is displayed in 
Fig 6. When the upper vortex has a clockwise circulation (F^, = — 1), the lift 



of the airfoil, produced by interaction, is negative ioi y 1 + h + h < — 2F^ 
or h < |, and positive for /i > |. Hence there is a distance at which the 
fiat plate has no lift gain from the upper vortex. When the upper vortex has 
anticlockwise circulation (F^, = 1), the lift is always positive. 



cl 
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Figure 6: Lift coefficient of the flat plate, which is interaction with a bound 
vortex at h above the middle point of the plate. 
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4 Summary 

Started from a momentum balance analysis, and proceeded with inter ex- 
change between the singularity velocity and induced flow velocity, we have 
obtained force formulas in both singularity form ( see (1121) for method of 
singularity velocity, and ( iT9l) or ( l24l) for single body or multibody method of 
induced velocity) and integral form (see f l25|) for method of singularity ve- 
locity and f l26l) for multibody method of induced velocity), for which bound 
singularities, multiple free singularities and multiple bodies can be consid- 
ered. The influence of the adjacent bodies on the actual body can be simply 
treated as influence of singularities representing the adjacent bodies, in a 
similar way as free singularities. Moreover, the influence on the force by 
vortex production is treated in a simple and explicit way (see ( 123|) and (!28|) . 
or Remark 2.2 in section 2.6). 

The present work is new for four reasons. First, it covers the work of 
Wu, Yang &Young (2012) as a special case, and includes in addition the 
effect of bound vortices and vortex production (see Remark 2.3 in section 
2.6). Second, the way to obtain the force formulas is based on the interac- 
tion of various singularities, and thus is useful for explicitly interpreting the 
influence due to various resources. For instance, we have shown that the 
interaction between free singularities do not contribute to forces, while the 
induced velocity effect is due to interaction between free singularities and 
inner singularities (see section 2.3). Third, the present result includes the 
situation where the discrete singularities are replaced by a distribution of 
vortices and source doublets, thus we have an integral approach which, with- 
out the use of an auxihary function, gives individual forces for each body 
in the case of multiple bodies (see section 2.5). Last, the present study ap- 
pears to provide a bridge between the singularity velocity approach, induced 
velocity approach and some integral approaches. 

The present results appear to be useful for deriving analytical force for- 
mulas even when multibody is considered. The validation study presented in 
section 3 demonstrates that analytical force formulas can be obtained even 
for rather complex problems. The example for an airfoil on top of which 
there is another airfoil (actually represented by a bound vortex) shows that 
it is possible to use the present method to derive analytical interaction forces 
for multiple airfoils, and this will be for a future study. 

The present results are actually restricted to two dimensional inviscid 
flow. For viscous flow it appears that the viscous effect can be treated sepa- 
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rately, by simply adding an additional boundary integral (Howe 1995). More- 
over, we did not consider the influence of body deformation, for which well es- 
tablished theories exist (Landweber&Miloh 1980, Howe 1995, Kanso&Oskouei 
2008). Finally, the case with body acceleration and rotation has not been 
tested in this paper. These problems need further studies. 

Similar analysis for axisymmetric flow has been carried out for insect 
flights (Wang&Wu 2010, Wang&Wu,2012), for which the influence of free 
vortices (in the form of vortex rings, or axisymmetric Karman vortex street) 
and insect body on the hft produced by the flapping wing is related to the 
induced velocity. In these studies the influence of the speeds of image vortices 
has been neglected, hence requiring further studies. 

A Momentum change inside a fixed body 

Now will show that 

when the body A is considered fixed. 

For each vortex or source i inside the body, we define an infinitesimal fixed 
circle of radius r^ = a/(x — Xj)^ + iy — ViY = ^ whose center instantaneously 
coincides with the vortex. For convenience we use here the subscript x, y and 
t denote the partial differentials with respect to x,y and t. Let E be the region 
bounded by the contour of the body dA and the circumferences of the circles 
Ti = e. Below are some additional relations. 

Decompose La defined in f lBU]) as 

La= P-f. I I vdydx + P^-ji I / vdydx 

Since v = —ipx and since ip is analytical in S, we use the divergence theorem 
to write 

/ / vdydx = ~ ipxdydx = — f tpdy + y^ (p ipdy 

Hence 

La = -jjHy + j:j^fHy + pj:jJl^^vdydx 



dA '•'" n=s 
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Consider the last term on the right hand side. Across an element of 
length Tiddi at (rj, 9i) on the circumference of the circle Tj < e, the loss of 
momentum due to the movement of the vortex (?) is 

drrii = -v^'^ {xi^t cos 6i + yi^t sin 6i) ndOi 
Remark that v^''^ = ^^ and 

-T, j I vdydx = ^ (^'^i = ^'^l7~ cos 9i{xi^t cos 9i + yi^t sin 9i) d9i 

•^ •^^i<^ i,in -^^ i,in "^^ 

Hence 

dec . . s-^ r. 



i,in %,in 

If sources can be similarly analyzed. When both vortices and sources are 
present we may write 

i,in >' >' I ^t j^j„ \ / 

Now consider the second term on the right hand side. First consider vortices. 

Let (x, y) be a fixed point on the circle. Differentiating x — Xi = ri cos 9i and 
y — yi = Ti sin 9i with respect to time we obtain 

ri9i^t = Xi^t sin 9i - y-i^t cos 9i, Ti^t = -Xi,t cos 9i - yi^t sin 9i 

Remark also that j^"" cos^ 9^^ d9^^ = n and /^^^ cos^« sin^Wrf^W = 0. Thus 
and similarly we have 

"p /'27r "p 

-^^dy = --^ / ri^tCos9id9i = -^x^ 
Ztt Jq I 

ri=e 



Hence 



iAn _ _ iAn 
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When sources are included, the analysis is similar and we may write 



-1 -iV) -1 -in ^ 



mi 



r, -^i,* r) yi,t 



In summary we have proved 



dA 

Since here the body is assumed stationary so that ip is & constant along the 
body, and thus § ipdy = 0. This means La = 0. 

dA 

Similarly we may prove Da = 0. 



B Additional expression for Karman vortex 
street 

To prove (l34l) using the present force formula flT2!) . we need a relation be- 
tween the velocities of the inner and outer vortices. For this purpose will 
define a large contour dC enclosing the body and a part of outer vortices. 
Moreover, for each outer vortex inside dC , we define an infinitesimal fixed 
circle of radius Vi = e whose center instantaneously coincides with the vortex. 
Now consider the fluid region F enclosed by dC, dA and A, where A denote 
the perimeters of all the fixed circles Vi = e. We derive some integrals along 
the contours dF- 

Since and ip are analytical in F, we may use the divergence theorem 
and the identity (j)y = —ipx , 4>x = i'y to write 

{(f)dx — ipdy) = ~ (01/ + V'x) dxdy = 

dF 

{(jydy + '(/'dx) = {4>x — "ipy ) dxdy = 

dF 



Hence 



^ {ipdy — (f)dx) — § {(pdx — ipdy) = § {ipdy — cpdx) 

dA A dC 

§ {ipdx + (pdy) + § {4>dy + ipdx) = § {4>dy + ipdx) 
dA A dc 
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(62) 



Similarly, for each vortex inside dA, we define an infinitesimal fixed circle 
of radius Vi = e whose center instantaneously coincides with the vortex. Now 
consider the region S enclosed by dA and B, where B denote the perimeters 
of all the fixed circles Ti = e inside dA. Since and if) are analytical inside 
the region E, we can apply the divergence theorem to § {^dy — (pdx) to write 

dA 

(b {ipdy — (f)dx) — (b {i^dy — (f)dx) = {(f'y + ^x) dxdy = 

dA e 

(b {ipdx + (t)dy) — (b {ijjdx + (j)dy) = (0x ~ "^y) dxdy = 

dA e 

Hence 

ii.dy - Hx) = / m - Hx) . / (Hx + ^dv) = / (i,dx + Hy) 
dA e dA e 



With the above relations we may rewrite ( l62l) as 

§ {ipdy — (pdx) = § {ipdy — (pdx) — § {(pdx — ipdy) 
dc e A 

§ {4>dy + ipdx) = § {ipdx + 4>dy) + § {4>dy + ipdx) 
dc A 

As for ( 16T|) in Appendix A, we may similarly show that 
§ (Ptdx = - § ^tdy = -^^ 

f My= § ^tdx = -Yj 

Inserting these expressions into (1631) . we obtain 



(63) 



dC 
dC 



(64) 



where Yli in ^^ ^^^^ vortices inside A and Yli f is over vortices inside F- 
The force formula ( !T2l) is split here as 
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i,in i,ou 



^ - -pEr^t-pEr 



which, when using flMj) to replace the first terms on the right hand side, 
yields 

L = pi § iijdy - <pdx) - pZr^^ 

dC i,D /^c;^| 

D = pfJ{(Pdy + ^|Jdx)+pET.^ ^ ^ 

aC i,D 

Here I? denotes the region outside of the contour C. 

Now, we introduce a downstream boundary x = x^ and assume that this 
is the contour dC. Then with (|H5]) we may write 



"""sL^'^^"..?/-^ 



^-^d i,Xi>xa 



Through defining ijj{x,y) = J2^j^i^ ~ ^j^V ~^ Vj) (where S is the Dirac 

j 
function), we may write 



E dyi d f 



i,Xi>xa Jx>xc 



With the identity yu = V ■ {yv, —yu) + m, and the divergence theorem so 
that /^^^, V ■ {yv, —yu) dxdy = — f^^^ yvdy, we further have 

^ Ti— ^ = p— y . [yv ^ -yu) dxdy + p— j udxdy 

i,Xi>xa •^^>^d "''' Jx>xa 

= —p-f / yudy + p— / udxdy 

"t- J X=X^ '^^ J X>Xa 

= P-T. I ^^y + P'i; / udxdy 

"''' J X=Xci "''^ Jx>Xd 

where we have used f yvdy = f yctiydy = — f (hdy. Hence 
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d f , , d f 

D = 2p— / (pay + p— / udxdy 

"t- J X=Xa "'^ J X>Xa 

Using the integral form of the y momentum equation for x > Xd, we have 
p-i: / udxdy = / (pw^ +p- pul^ - Poo) dy 

^'^ Jx>X^ J X = Xc 

and with the BernouUi equation p = —^{u^ + v^) — (pt + Poo + I^L^ "^^ have 
^ ^ Pjt ^^y+nP ("^ -v^ - M^) dy 

11- Jx=x,t ^ J x=x^ 



Since the above analysis is frame independent, we may choose a frame 
attached to the vortex street and therefore 0j = on the line x = Xd, 
which is assumed to intersect the vortex street in its uniform region. Hence 
D = ^p f^^^ {v? — v^ — M^) dy. Now we decompose (m, v) qs u = Mqo + Wfc, 
^ = ■^00 + ^^ = ^fc, where (m^, f jt) is the induced velocity by the vortex street, 
then 

^ = \p [ H-4)dy + Dr 

^ Jx=xa 

where Dr = puoo f _ Ukdy. It is obvious that Dr = since the contribution 
to Uk by any vortex is antisymmetric about the y position of this vortex. 
Thus we proved ( 1M|) in section 3.2. 

C Additional expression for the Wagner prob- 
lem 

Remark that for the Wagner problem 

1 r 2 r 

Aq = -a -— / v^{^)dp, An,n>i = —rj- / v^{C) cos (ra/3) ci/3 

where v^{^) is induced by the vortex sheet and is given by 

"^^ 27rJo CA + s-^ 27rJo f(cos/3 + l) + s 
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Thus 



^0 = -« + ^^ lo^' Jo f^(Jg+i)+, fe(g, t)ds 

^2 - 2^2yo„ Jo Jo Jo i2L(cos/3+l)+s'^^^'^^^^ 



Hence 



1. 1 f^°°' r 2 (cos^ /3 - COS /3) d/3 ^ , ,, 

It is straightforward to show that 

' ^0 + lA, = -« - ^^ /o^-* (1 - ^) M^, t)rf. 



'^"' 2 ttVocca Jo Vs 



by using 



rn 

Jo 3a 



d/3 



2-(cos/3+l)+s v^\/cA+i 

cos/3-1 ^ _2 2(eA+5) 

^(cos/3+l)+s CA ca(^(cos/3+1)+s) 



2 
cos^/3-cos/3 _ 2COS/3 _ 4(ca+s) , 2(cA+2g)(cA+s) 

^(cos/3+l)+s ~ CA ~~^ c^(fA(cos/3+l)+s) 
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